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SUMMARY 

This  paper  describes  changes  to  the  computer  program  "SUPER", 
which  correct  errors  identified  in  the  previous  version  of  the  program,  and 
enable  the  calculation  of  generalised  forces.  Tables  of  generalised  forces, 
calculated  using  "SUPER"  and  other  programs,  are  presented  for  a  standard 
AGARD  configuration. 
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1.  INTRODUCTION 


Recently,  a  computer  program  (SUPER)  has  been  developed  for  the  calcu¬ 
lation  of  linearised  steady  and  oscillatory  supersonic  flow  over  a  configuration 
consisting  of  a  coplanar  wing  and  tail  in  association  with  a  vertical  fin  [1],  This 
program  implements  an  explicit  finite  difference  scheme  based  on  [2],  extending 
and  correcting  work  done  previously  at  ARL  ([3]).  The  primary  motivation  for 
this  work  is  to  enable  the  calculation  of  unsteady  aerodynamic  derivatives,  i.e. 
generalised  forces,  for  insertion  into  aeroelastic  mathematical  models. 

In  SUPER,  the  velocity  potential  is  calculated  throughout  the  flow  field,  and 
both  the  potential  and  pressure  are  calculated  on  the  lifting  surfaces.  However, 
the  program  stops  short  of  calculating  generalised  forces  from  these  quantities. 

The  purpose  of  this  paper  is  to  list  and  correct  a  number  of  errors  identified  in 
the  last  version  of  SUPER  (i.e.  in  [1]),  and  to  describe  the  modifications  to  SUPER 
which  enable  the  calculation  of  the  required  generalised  forces.  Generalised  forces 
calculated  using  SUPER  are  compared  with  the  results  of  other  solution  schemes 
for  a  number  of  standard  configurations. 


2.  REVISED  DESCRIPTION  OF  PROGRAM 


Only  significant  changes  to  SUPER  will  be  described  in  this  section,  e.g. 
changes  which  alter  the  input  to  or  output  from  the  program,  or  which  correct 
errors  in  the  program. 

The  old  “wing  functions”,  “tail  functions”  and  “fin  functions”  parameters 
IPAR,  JPAR  and  KPAR  respectively  have  been  replaced  by  a  single  “forcing  mode” 
number  IFM.  The  last  version  of  SUPER  could  calculate  the  potential  and  pressure 
for  a  regular  sequence  of  Mach  numbers  and  reduced  frequencies,  and  only  for  a 
single  “forcing  mode” .  Modifications  have  now  been  made  to  allow  the  program 
to  handle  a  list  of  Mach  numbers,  reduced  frequencies  and  “forcing  modes”.  The 
“forcing  modes”  have  been  further  subdivided  into  symmetric  and  antisymmetric 
groups,  so  that  the  generalised  forces  for  the  two  groups  can  be  calculated  and 
displayed  separately. 

Calculation  of  the  generalised  forces  has  naturally  led  to  additions  and  modifi¬ 
cations  to  the  code.  There  are  now  three  new  subroutines  WINGFORCES,  TAIL- 
FORCES  and  FINFORCES,  in  which  contributions  to  the  generalised  forces  are 
calculated.  Output  from  these  subroutines  is  collated  in  a  new  array  GENFOR 
which  contains  the  generalised  forces  from  each  mode/Mach  number/frequency 
run.  The  symmetric  and  antisymmetric  results  are  later  extracted  from  GENFOR 
for  tabulation. 

The  subroutines  SOLVEI  and  SOLVE2  have  been  changed  to  prevent  discrep¬ 
ancies  at  the  leading  and  trailing  edges  between  these  routines  and  the  output  rou¬ 
tines  OPUT1  and  OPUT2.  These  changes  entail  working  with  the  indices  (i,j,k) 
of  each  point,  as  in  OPUT1  and  OPUT2,  instead  of  the  cartesian  co-ordinates 
(x,y, a),  when  implementing  boundary  conditions. 
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Two  errors  were  discovered  and  corrected  in  SOLVE1  and  SOLVE2.  In  the 
program  listing  on  page  84  of  [lj  the  statement 

Z1=FL0AT(K-1)/DD 

appears.  The  correct  line  is 

Zl=FLOAT(K-i)/DD/2.0  . 

Tests  using  the  standard  AGARD  planform  (see  [5])  revealed  this  error,  which 
was  not  detected  in  [1]  because  the  code  was  tested  only  for  a  rectangular  fin 
(for  which  the  leading  and  trailing  edge  positions  are  independent  of  the  vertical 
coordinate  Zl). 

The  second  error  involved  setting  the  contribution  from  the  derivative  bound¬ 
ary  condition  points  to  zero  when  one  or  both  of  these  points  lay  off  the  lifting 
surface  (in  a  wake  or  diaphragm  region)  immediately  ahead  of  a  solution  point 
lying  on  the  surface.  These  normal  derivatives  are  in  general  non-zero  and  un¬ 
known  in  advance.  A  quick  and  consistent  solution  to  this  problem  was  found  to 
be  evaluation  of  the  derivatives  as  if  the  points  were  actually  on  the  lifting  surface. 

With  the  above  modifications,  input  to  the  program  has  changed  correspond¬ 
ingly.  A  new  value  of  the  output  parameter  IOUT  has  been  added  (IOUT=4) 
to  enable  the  generalised  forces  to  be  printed  without  the  accompanying  pressure 
and/or  potential  output.  As  mentioned  above,  the  variables  IPAR,  JPAR  and 
KPAR  have  been  replaced  by  IFM. 

Multiple  “forcing  modes”,  Mach  numbers  and  reduced  frequencies  are  handled 
by  writing  the  number  of  each  quantity  on  one  line,  and  the  list  of  values  in  free 
format  on  the  line,  or  lines  below. 

The  “displacement  modes”  to  be  considered  are  separated  according  to  their 
spanwise  symmetry.  The  spanwise  symmetric  modes,  of  which  there  are  NSM,  are 
numbered  first,  followed  by  the  NAM  antisymmetric  modes,  which  are  numbered 
NSM+1  to  NSM+NAM.  Each  “forcing  mode”  must  belong  to  one  of  these  two 
subsets,  so  that  its  symmetry  is  clear  to  the  program. 

Finally,  the  streamwise  grid  spacing  has  been  replaced  by  the  spanwise  grid 
spacing  as  input  to  the  program.  (One  can  easily  be  calculated  from  the  other.) 

The  new  input  file  takes  the  form 

1  1  if  flow  unsteady,  2  if  flow  steady. 

2  Output:  1  press.,  2  pot.,  3  both,  4  gen.  forces 

1  Width  :  1  for  80  char,  wide  output,  2  for  132  char.. 

6  No.  of  dec.  places  in  output,  must  be  in  range  4-8. 

0.05  Spanwise  grid  step  size. 

1.0  Wing  semi  span  length. 

1.0  Tail  semi  span  length. 

0.0  Lowest  vertical  point  of  fin. 

1.2  Highest  vertical  point  of  fin. 

2  Number  of  Mach  no . s 

1.2  3.0  Mach  numbers 

3  Number  of  Red.  freq. 

1.5  0.5  1.0  Reduced  frequencies 
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0  No.  of  symm.  displ.  modes  (labelled  1...NSM) 

4  No.  of  antisymm.  displ.  modes  (NSM+1 . . .NSM+NAM) 

2  Number  of  forcing  modes  (NFM) 

1  3  Forcing  modes. 

The  separation  of  the  modes  into  symmetric  and  antisymmetric  groups  has 
also  led  to  necessary  changes  in  the  mode  shape  subroutines.  For  example,  corre¬ 
sponding  to  the  input  data  file  above,  these  subroutines  take  the  form  shown  in 
Appendix  1. 

The  first  mode  shape  in  each  subroutine,  corresponding  to  label  1,  refers  to 
the  shape  of  the  wing  for  the  steady  flow  calculation.  Subsequent  labels  refer 
to  the  unsteady  modes,  e.g.  label  2  corresponds  to  unsteady  mode  1.  Looking 
“horizontally”  across  the  mode  shape  subroutines,  it  can  be  seen  that  mode  1  is 
purely  a  wing  mode  (twist)  since  there  are  zeros  in  the  equivalent  positions  in 
the  tail  and  fin  subroutines.  Similarly  mode  2  is  a  pure  tail  mode  (pitch),  and 
modes  3  and  4  are  purely  fin  modes  (bending  and  twist  respectively).  No  doubt 
each  subroutine  could  be  shortened  by  eliminating  all  but  one  of  the  zero  modes; 
however  this  would  contribute  nothing  to  the  clarity  of  the  code. 

3.  CALCULATION  OF  GENERALISED  FORCES 


SUPER  computes  the  potential  and  pressure  on  each  lifting  surface.  There 
are  two  approaches  to  the  calculation  of  generalised  forces  from  these  quantities. 
The  first  of  these  involves  direct  integration  of  pressure  multiplied  by  mode  shape 
over  the  whole  surface  (which  to  fix  ideas  we  will  take  here  to  be  the  wing). 
The  generalised  force  Qy  in  the  »'th  mode  of  oscillation  resulting  from  harmonic 
oscillation  in  the  j th  mode  is  given  by 

t  lie(y) 

Qa  =  J  J  ACPJ-(i, y)fi (x, y)dxdy.  (2.1) 

-t 

where  l  is  the  wing  semi-span, 

ACPJ(x,y)  =  -2(A<t>JX(x,y)  +  iuA<t>3(x,y)), 

is  the  difference  in  pressure  coefficient  between  the  lower  and  upper  surfaces  of 
the  wing  resulting  from  oscillation  in  mode  j  (the  jth  “forcing  mode”),  fi{x,y)  is 
the  tth  “displacement  mode”  and  w  is  the  reduced  frequency.  The  above  formula 
is  equivalent  to  that  given  in  [4] .  Similar  formulae  apply  for  the  tail  and  fin. 

This  straightforward  method  involves  use  of  both  the  reduced  potential  tj>,  (ac¬ 
tually  4> o  in  [  1]) ,  and  its  streamwise  derivative  <j>x.  There  are,  however,  significant 
problems  associated  with  the  use  of  4>x,  and  these  are  described  below. 

There  is  a  fundamental  difference  between  the  behaviour  of  the  pressure  dis¬ 
tribution  at  a  subsonic  leading  edge,  where  it  has  a  square-root  singularity  and  at 
a  supersonic  leading  edge,  where  it  remains  bounded.  For  the  generalised  force  re¬ 
sults  to  be  accurate,  the  integration  scheme  must  determine  whether  a  singularity 
exists  and  model  it  correctly. 


This  is  particularly  difficult  because,  in  general,  the  streamwise  potential  con¬ 
tains  jagged  oscillations.  These  oscillations  arise  because  the  geometry  of  the 
leading  edge  of  the  wing  is  not  modelled  exactly,  appearing  to  the  program  to 
be  changing  in  a  steplike  manner.  A  pair  of  Mach  lines  emanate  from  each  of 
these  step  changes,  and  sharp  changes  occur  in  the  streamwise  potential  as  each 
line  is  crossed.  On  differentiation  these  oscillations  are  magnified,  swamping  the 
true  behaviour  of  the  pressure  at  the  leading  edge,  and  making  it  nearly  impos¬ 
sible  to  integrate  accurately.  As  a  result,  the  direct  integration  method  can  only 
confidently  be  applied  to  supersonic  leading  edges. 

The  effect  of  the  jagged  leading  edge  has  been  noted  previously  (e.g.  in  [12]) 
in  conjunction  with  the  Mach  box  method,  and  the  solution  has  been  to  bypass 
the  need  for  the  pressure  by  carrying  out  a  preliminary  integration  by  parts.  This 
gives  the  generalised  force  integrals  in  terms  of  the  potential  only.  To  be  specific, 
it  leads  to  the  formula 

Qa=  73  /  A<M*»v)/<(*.v)  dv 
1  it  L  J  *..(») 

/  *!•(») 

~jaj  J  A^j(x>  S') 

-i  ii.(v) 

The  advantage  here  is  that  since  <j>j  is  always  bounded,  it  requires  no  special 
treatment  (in  the  surface  integral  at  least).  One  possible  disadvantage  of  this 
method  is  that  it  involves  line  integration  at  both  the  leading  and  trailing  edges  and 
these  could  conceivably  be  sensitive  to  the  streamwise  oscillations  in  <f>.  However, 
experience  shows  that  the  effect  of  the  oscillations  is  reasonably  small  in  the  line 
integration. 

As  before,  the  main  problem  lies  with  subsonic  leading  edges,  i.e.  regions  in 
which  <j)j  changes  rapidly.  It  was  found  that  the  resolution  at  the  leading  edge 
was  not  always  adequate  to  capture  the  full  effect  of  the  singularity  in  slope.  This 
was  remedied  by  evaluating  in  the  line  integral  at  a  position  one  grid  point 
upstream  of  the  leading  edge.  The  philosophy  here  is  that  the  potential  changes 
relatively  slowly  ahead  of  the  lifting  surface,  so  the  error  in  estimating  a  rapid 
change  in  <j>}-  near  the  leading  edge  should  be  small. 

For  the  wing  and  fin,  the  above  solution  can  always  be  applied  because  the 
potential  is  identically  zero  in  the  plane  of,  and  ahead  of,  the  wing  or  fin.  For  a  tail 
immersed  in  the  wake,  this  solution  should  apply  as  long  as  the  reduced  frequency, 
which  controls  the  frequency  of  oscillation  of  4>  in  the  wake,  is  not  too  large. 

In  the  latest  version  of  SUPER  we  implement  the  second  method  described 
above  for  the  calculation  of  generalised  forces.  All  integrals  are  performed  using 
the  trapezoidal  rule  in  one  and  two  dimensions.  The  results  presented  in  Section 
4  show  reasonable  convergence  properties  with  decreasing  grid  size. 


4.  RESULTS  AND  DISCUSSION 


Comparison  has  been  made  between  the  output  from  SUPER  and  that  from 
a  number  of  other  methods,  namely 

a.  the  Mach  box  method,  [6j  and  [7], 

b.  the  potential  gradient  method,  [8]  and  [9]  (improved), 

c.  the  doublet  point  method,  [10],  and 

d.  the  harmonic  gradient  method,  [11], 

for  a  number  of  well  known  interference  configurations.  In  addition,  results  are 
also  given  from  [2]  in  which  the  present  finite  difference  scheme  first  appeared. 

It  should  be  noted  here  that  [4]  provides  an  exhaustive  list  of  generalised  force 
results  for  a  number  of  wing-only  planforms.  However,  as  agreement  with  those 
results  is  quite  good,  only  those  results  for  the  more  difficult  cases  of  coplanar 
wing  and  tail,  with  and  without  vertical  fin  interference  are  presented  here.  As 
mentioned  above,  the  test  cases  used  are  largely  defined  in  the  AGARD  supplement 
[5]  (from  which  Figure  1  is  reproduced),  and  many  of  the  results  quoted  here 
appear  in  that  reference.  Forty  spanwise  grid  points  were  tried  as  in  [2],  and  gave 
reasonable  results.  However,  it  was  decided  to  generate  the  present  results  with 
sixty  spanwise  points,  as  this  led  (naturally)  to  better  agreement,  especially  at 
M  =  3.0. 

Tables  1  to  3b  show  generalised  forces  calculated  for  the  coplanar  wing  and 
tail  configuration  shown  in  Figure  1  at  Mach  numbers  1.2,  1.356,  and  3.0,  and 
reduced  frequencies  0.0  and  1.5.  The  results  of  the  present  method  (which  here 
and  below  we  will  refer  to  as  [l])  are  compared  with  those  of  [2],  [6],  [7],  [8]  and 
[10].  The  mode  shapes  here  are 


1.  Wing  Twist 

2.  Wing  Bending 

3.  Tail  Roll 
and  4.  Tail  Pitch 


z  =  y[x  —  2.25|y|  —  0.85), 

«  =  y|y|, 

2  =  y, 

z  =  (x~  3.35)sgn(y). 


There  is  good  agreement  here  between  the  various  codes,  except  that  the  present 
method  tends  to  underestimate  the  magnitude  of  the  generalised  forces  by  a  small 
amount. 

Table  4  shows  generalised  forces  obtained  from  the  present  program  for  a 
coplanar  wing,  tail  and  vertical  fin  at  Mach  number  3.0.  Results  from  the  present 
method  for  this  configuration  are  compared  to  results  from  [8]  and  [9]  for  the 


following  modes: 

1. 

Wing  Twist 

z  =  y(x  -  2.25|y|  —  0.85), 

2. 

Tail  Pitch 

2  =  (x  -  3.35)  sgn  (y), 

3. 

Fin  Bending 

y  =  21 2 3, 

and 

4. 

Fin  Twist 

y  =  z( i  -  0.875z  —  3.0). 

The  results  of  [9]  appear  (as  expected)  to  be  more  reliable  than  those  of 
[8].  Certainly,  there  is  better  agreement  in  the  orders  of  magnitude  of  the  terms 
between  [1]  and  [9].  Notice,  however,  that  there  are  some  significant  differences  in 
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sign  and  magnitude  for  the  modes  involving  perpendicular  surfaces  e.g.  Q3l  (fin 
and  wing)  and  Q24  (tail  and  fin),  while  the  generalised  forces  for  modes  involving 
the  same  or  coplanar  surfaces  all  agree.  This  suggests  a  possible  error  in  the 
solution  method  of  either  [1]  or  [9]  (as  opposed  to  the  post-processing  to  compute 
the  generalised  forces).  The  solution  method  in  [l]  has  been  checked  in  a  number 
of  ways  to  ensure  consistency  between  the  y  and  z  directions,  and  for  the  correct 
sign  of  the  generalised  forces  for  some  simple  modes. 

Results  in  Tables  5a  and  5b  show  similar  behaviour  to  that  shown  in  Table 
4.  Here  comparison  is  made  between  the  present  method  and  the  results  of  [8], 
[9]  and  [ll]  for  the  T-tail  configuration  shown  in  Figure  1.  The  present  program 
is  set  up  so  that  the  horizontal  tail  is  in  the  plane  z  =  0  instead  of  z  =  1.2.  The 
mode  shapes  are  adjusted  accordingly  to  become 

1.  Fin  Bending  y  =  (2  +  1.2)2, 

2.  Fin  Twist  y  =  (z  +  1.2)(x  —  0.875(z  +  1.2)  —  3.0), 

and  3.  Tail  Roll  z  =  y. 

Again,  there  are  significant  differences  between  the  codes  for  modes  involving 
perpendicular  surfaces,  while  the  modes  involving  coplanar  surfaces  are  in  good 
agreement. 


5.  CONCLUSION 


This  paper  serves  two  purposes.  The  first  is  to  describe  changes  to  the  pro¬ 
gram  SUPER  presented  in  [l]  resulting  from  the  correction  of  programming  errors 
and  minor  inconsistencies,  and  from  alterations  enabling  the  calculation  of  gen¬ 
eralised  forces.  The  second  is  to  compare  the  generalised  force  results  with  those 
calculated  by  different  methods  for  some  standard  test  cases. 

It  is  found  that,  in  general,  the  results  from  the  present  method  agree  quite 
well  with  comparable  codes,  except  that  the  present  method  tends  to  underesti¬ 
mate  the  magnitude  of  the  forces  by  a  small  amount.  However,  there  are  some 
differences  for  off-diagonal  modes  involving  perpendicular  surfaces  (i.e.  tail  and 
fin,  wing  and  fin)  which  cannot  fully  be  explained  at  present.  In  those  cases, 
though,  there  is  some  scatter  in  the  results  of  all  methods.  Careful  testing  of  the 
present  code  has  produced  no  evidence  of  internal  inconsistency  or  error. 

Future  work  with  SUPER  will  involve  computing  generalised  forces  for  real¬ 
istic  geometries  such  as  the  Fill  (wing,  laii  and  vertical  fin,  with  the  fuselage 
represented  by  horizontal  and  vertical  surfaces). 
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TABLE  2 


Wing  and  Horizontal  Tail  Interference 
at  M  =  1.356  and  k  =  0.0, 1.5 


■ 

?r 

II 

O 

o 

k  =  1.5 

ns 

Ref.  6 

Ref.  7 

Ref.  1 

Ref.  6 

Ref.  7 

Ref.  1 

Q'a 

i,i 

-0.0541 

-0.0587 

-0.0484 

0.0041 

-0.0552 

2,1 

0.2910 

0.2859 

0.4637 

0.4986 

0.4226 

3,1 

-0.5073 

-0.4986 

-0.3395 

-0.4503 

-0.3258 

4,1 

-0.3218 

-0.3024 

-0.1742 

-0.2507 

-0.1671 

1,2 

-0.1501 

-0.1523 

-0.1361 

2,2 

-0.2831 

-0.2694 

-0.2623 

3,2 

-0.4369 

-0.4341 

-0.4083 

4,2 

-0.2882 

-0.2941 

-0.2693 

1,1 

0.2766 

0.2992 

0.2461 

2,1 

0.3933 

0.4168 

0.3705 

3,1 

0.3261 

0.3115 

0.3021 

4,1 

0.2083 

0.2083 

0.1925 

1,2 

0.0120 

0.0299 

0.0063 

2,2 

0.2329 

0.2912 

0.2667 

3,2 

0.0160 

-0.0144 

0.0116 

4,2 

0.0179 

-0.0041 

0.0132 
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TABLE  3a 

Wing  and  Horizontal  Tail  Interference 
at  M  =  3.0  and  k  =  0.0 


TABLE  4 


TABLE  5a 

Tail  and  Fin  Interference 
M  =  1.6,  Z  =  1.2  and  ib  =  0.0 


TABLE  5b 


Tail  and  Fin  Interference 
at  M  =  1.6,  Z  =  1.2  and  k  =  1.5 


■ 

m 

k  = 

1.5 

Ref.  8 

Ref.  9 

Ref.  11 

Ref.  1 

mm 

0.0093 

mm 

0.0298 

HI 

pi 

-0.0553 

-0.0509 

Hi 

0.0560 

0.0473 

m 

-0.3310 

MU 

mm 

0.6926 

0.6864 

0.6879 

1 

0.1158 

0.1250 

0.1159 

|| 

-0.1359 

-0.0458 

-0.0492 

K 

MEM 

0.0450 

0.0676 

-0.0322 

IH 

El 

0.1795 

0.0277 

-0.0132 

0.0192 

0.0622 

0.0266 

Q", 

mm 

0.5200 

0.5457 

0.5466 

0.5444 

mm 

0.0449 

0.0462 

0.0620 

0.0462 

B9 

0.1209 

-0.0627 

0.0471 

-0.0373 

mm 

0.1295 

181 

0.1518 

0.1497 

mm 

■SEffli 

-0.0071 

-0.0952 

0  2091 

1,3 

0.0527 

0.0341 

-0.0268 

2,3 

0.0166 

0.0157 

0.0074 

-0.0076 

3,3 

0.4457 

0.4783 

0.4309 

0.4188 
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APPENDIX  1  -  LISTING  OF  MODE  SHAPE  SUBROUTINES 


C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

REAL  FUNCTION  WL(Y) 

WL=2.75*Y 
RETURN 
ENTRY  ViT(Y) 

WT=2 .25+1 .45*Y 
RETURN 
ENTRY  TL(Y) 

TL=2 .7+1 .2*Y 
RETURN 
ENTRY  TT(Y) 

TT=4.0+0.26*Y 
RETURN 
ENTRY  FL(Z) 

FL=2.6+Z 
RETURN 
ENTRY  FT(Z) 

FT=4.2+Z/3.0 

RETURN 

END 

C 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE  WING(X.Y.IP,H1,H2,IER) 

IER=0 

GO  TO  (1.2. 3. 4, 6)  IP 

IER=1 

RETURN 

1  CONTINUE 
Hl=-1 .0 
H2=-l  .0 
RETURN 

2  CONTINUE 

H1=Y* (X-2 . 25*ABS (Y) -0.85) 

H2=Y 

RETURN 

3  CONTINUE 
H1=0.0 
H2=0.0 
RETURN 

4  CONTINUE 
Hl-0.0 
H2=0.0 
RETURN 

6  CONTINUE 

H1=0.0 
H2=0.0 
RETURN 
END 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 
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c 

cccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccccc 

c 

SUBROUTINE  TAIL(X, Y, IP ,H1 ,H2, IER) 

IER=0 

GO  TO  (1.2. 3, 4, 5)  IP 

IER=1 

RETURN 

1  CONTINUE 
Hl=-1 .0 
H2=-l .0 
RETURN 

2  CONTINUE 
H1=0.0 
H2=0 . 0 
RETURN 

3  CONTINUE 
H1=(X-3.36)*SIGN(1.0.Y) 

H2=SIGN(1 .0. Y) 

RETURN 

4  CONTINUE 
H1=0.0 
H2=0.0 
RETURN 

6  CONTINUE 

H1=0.0 
H2=0 . 0 
RETURN 
END 
C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 

c 

SUBROUTINE  FIN(X .Z .IP. HI ,H2 ,IER) 

IER=0 

GO  TO  (1.2, 3, 4. 6)  IP 

IER=1 

RETURN 

1  CONTINUE 
Hl=-1 .0 
H2=-1.0 
RETURN 

2  CONTINUE 
H1=0.0 
H2=0.0 
RETURN 

3  CONTINUE 
H1=0.0 
H2=0.0 
RETURN 

4  CONTINUE 
H1=Z*Z 
H2=0.0 
RETURN 

5  CONTINUE 
Hl=Z*(X-0. 875*Z-3 .0) 

H2=Z 

RETURN 

END 

C 

CCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCCC 


17 


DISTRIBUTION 


AUSTRALIA 

Department  of  Defence 

Defence  Central 

Chief  Defence  Scientist 

FAS  Science  Corporate  Management  (shared  copy) 

FAS  Science  Policy  (shared  copy) 

Director,  Departmental  Publications 

Counsellor,  Defence  Science,  London  (Doc  Data  Sheet  Only) 
Counsellor,  Defence  Science,  Washington  (Doc  Data  Sheet  Only) 
S.A.  to  Thailand  MRD  (Doc  Data  Sheet  Only) 

S. A.  to  the  DRC  (Kuala  Lumpur)  (Doc  Data  Sheet  Only) 

OIC  TRS,  Defence  Central  Library 

Document  Exchange  Centre,  DISB  (18  copies) 

Joint  Intelligence  Organisation 

Librarian  H  Block,  Victoria  Barracks,  Melbourne 

Director  General  -  Army  Development  (NSO)  (4  copies) 

Defence  Industry  and  Materiel  Policy,  FAS 

Aeronautical  Research  Laboratory 
Director 
Library 

Chief  -  Aircraft  Structures 
Head  -  Flight  Mechanics  Branch 
Divisional  File  -  Aircraft  Structures 
Author:  I.H.  Grundy 
B.C.  Hoskin 
R.  Jones 
J.L.  Kepert 

T. G.  Ryall 
P.A.  Farrell 
B.  Emslie 
W.  Waldman 

Materials  Research  Laboratory 
Director/Library 

Defence  Science  &  Technology  Organisation  -  Salisbury 
Library 

Navy  Office 

Navy  Scientific  Adviser  (3  copies  Doc  Data  sheet  only) 

Army  Office 

Scientific  Adviser  -  Army  (Doc  Data  sheet  only) 

Air  Force  Office 

Air  Force  Scientific  Adviser 


Civil  Aviation  Authority 
Mr  C.  Torkington 


CSIRO 

Chemical  and  Wood  Technology  Division 

Universities  and  Colleges 
Adelaide 

Barr  Smith  Library 

Head,  Department  of  Mathematics 
ANU 

Head,  Department  of  Mathematics 
Canberra  CAE 

Head,  Department  of  Mathematics 
Chisholm  IT 

Head,  Department  of  Mathematics 
Flinders 

Head,  Department  of  Mathematics 
Footscray  IT 

Head,  Department  of  Mathematics 
Gippsland  IAE 

Head,  Department  of  Mathematics 
Griffith 

Head,  Department  of  Mathematics 
James  Cook 

Head,  Department  of  Mathematics 
Macquarie 

Head,  Department  of  Mathematics 
Melbourne 

Head,  Department  of  Mathematics 
Melbourne  IAE 

Head,  Department  of  Mathematics 
Monash 

Head,  Department  of  Mathematics 
Murdoch 

Head,  Department  of  Mathematics 
New  England 

Head,  Department  of  Mathematics 
NSW 

Library,  Australian  Defence  Force  Academy 

Head,  Department  of  Mathematics 
Newcastle 

Head,  Department  of  Mathematics 
Queensland 

Head,  Department  of  Mathematics 
QIT 

Head,  Department  of  Mathematics 
RMIT 

Library 

Dr  J.A.  Gear,  Department  of  Mathematics 

Head,  Department  of  Mathematics 
SAIT 

Head,  Department  of  Mathematics 


Swinburne  IT 

Head,  Department  of  Mathematics 
Sydney 

Head,  Department  of  Mathematics 
Tasmania 

Head,  Department  of  Mathematics 
Tas.  SIT 

Head,  Department  of  Mathematics 
Warmambool  IAE 

Head,  Department  of  Mathematics 
Western  Australia 

Head,  Department  of  Mathematics 
Woollongong 

Head,  Department  of  Mathematics 


SPARES  (10  copies) 
TOTAL  (88  copies) 


AL  149 


DEPARTMENT  OF  DEFENCE 


DOCUMENT  CONTROL  DATA 


PAGE  CLASSIFICATION 

UNCLASSIFIED 


PRIVACY  hARKING 


la.  AR  NUACER 

AR-005-591 


lb.  ESTABLISHMENT  NUMBER 

ARL-STRUC-TM-504 


2.  DOCTHQir  DATE 

APRIL  1989 


3.  TASK  NUMBER 

DST  87/023 


ADDITIONS  AND  CORRECTIONS  TO 
"SUPER"  A  PROGRAM  FOR  CALCULATING 
STEADY  AND  OSCILLATORY 
SUPERSONIC  FLOW  OVER  A  THIN  WING, 
TAIL  PLANE  AND  FIN 


8.  AITTHCR(S) 

I.H.  GRUNDY 


S.  SECURITY  CLASSIFICATION  6.  NO.  PAGES 

(PUCE  APPROPRIATE  CUSSIFICATION 
IN  BCK(S)  IE.  SECRET  (S),  COKF.(C)  18 

RESTRICTED  (R),  UHCLASSIFIED  (U)  ). 


□  HQ 


7.  HO.  REFS. 
12 


TITLE  ABSTRACT 


9.  DOWNGRADING /DELIMITING  INSTRUCTIONS 

Not  applicable. 


10.  CORPORATE  AUTHOR  AND  ADDRESS 

11.  CmCE/POSITION  RESPONSIBLE  FORt 

AERONAUTICAL  RESEARCH  LABORATORY 

DSTO 

SPONSOR 

P.0.  BOX  4331,  MELBOURNE  VIC  3001 

SECURITY 

DOWNGRADING 

CSTD 

APPROVAL 

12.  SBOCKMKT  DISTRIBUTION  (OT  THIS  DOCUMENT) 

Approved  for  public  release. 

cwnsms  bohries  outside  srxrm  umrKncm  should  ism  humid  through  asdis,  Duma  imfcmaxkm 

13a.  THIS  DOCUNBfT  MAT  «  i 

No  limitations. 


•  IN  CAIALOGCTB  1 


SWVICNS  AVAILABLE  TO. . 


13b.  CITATION  FOR  OTHER 
AMNOUMOMBfr)  HAT  BE 


PURPOSES  (IK.  CASUAL 


J  AS  FOR  13a. 


IS.  DRDA  SUBJECT 


00S1A 


Jt.  ABSTRACT 

This  paper  describes  changes  to  the  computer  program  *§UPER",  which 
correct  errors  Identified  In  the  previous  version  of  the  program, 
and  enable  the  calculation  of  generalised  forces.  Tables  of 
generalised  forces,  calculated  using  “SuPER^and  other  programs,  are 
presented  for  a  standard  AGARD  configuration. 


14.  DBC3UPT0RB 

Steady  flow 

Supersonic  flow 

Thin  wings 

Tall  assemblies 

Fins 

Oscillating  Flow 

Super  (computer  program) 

